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It has been recently claimed that the initial singularity might be avoided in the context of rain¬ 
bow cosmology, where one attempts to account for quantum-gravitational corrections through an 
effective-theory description based on an energy-dependent (“rainbow”) space-time metric. We here 
scrutinize this exciting hypothesis much more in depth than previous analyses. In particular, we 
take into account all requirements for singularity avoidance, while previously only a subset of these 
requirements had been considered. Moreover, we show that the implications of a rainbow metric 
for thermodynamics are more significant than previously appreciated. Through the analysis of two 
particularly meaningful examples of rainbow metrics we find that our concerns are not merely im¬ 
portant conceptually, but actually change in quantitatively significant manner the outcome of the 
analysis. Notably we only find examples where the singularity is not avoided, though one can have 
that in the regime where our semi-classical picture is still reliable the approach to the singularity 
is slowed down when compared to the standard classical scenario. We conclude that the study of 
rainbow metrics provides tantalizing hints of singularity avoidance but is inconclusive, since some 
key questions remain to be addressed just when the scale factor is very small, a regime which, as 
here argued, cannot be reliably described by an effective rainbow-metric picture. 

PACS numbers: 04.60.Be; 04.20.Dw; 98.80.Qc 


I. INTRODUCTION 

The avoidance of the initial singularity has been subject of intense debate in cosmology (see Ref. [I] and references 
therein). In particular, since we expect that general relativity (GR) should be no longer reliable under the extreme 
conditions of the primordial universe, where the curvature and/or its invariants diverge, quantum gravitational effects 
might provide a mechanism for avoiding the singularity. Examples of some attempts along this direction are given in 
Ref. 0. 

In order to investigate this possibility one of the most natural and manageable starting points is provided by the 
“rainbow gravity” picture. There one assumes that quantum gravity should admit a semi-classical regime in which 
quantum-gravitational corrections can be effectively described by introducing an energy-dependent spacetime metric, 
while otherwise keeping the framework rather conventional, with spacetime geometry coded in a smooth manifold. 
This was inspired by (and provides a link with) the intense effort devoted to the study [3HS] of the hypothesis of 
Planck-scale-modified dispersion relations (MDR): with an energy-dependent spacetime metric particles of different 
energies would experience different geometries of spacetime, possibly accounting for the expected fact that the short- 
distance quantum structure of spacetime would affect differently particles of different wavelengths (more strongly the 
ones of shorter wavelength). The notion of rainbow geometry was introduced in [7], as a perspective on the structure 
of “DSR-relativistic theories” IMD], and since then it was widely used for quantum-gravity-phenomenology analyses, 
see e.g. |1 ITtlSj . While at Erst the notion of rainbow metric had only a semi-heuristic interpretation, there are now 
formalisms that can consistently accommodate such a notion, such as the relative-locality framework JUKI?] and 
certain Finsler geometries using. 

We here propose to explore the possible role of rainbow metrics in singularity avoidance by following a full analysis 
centered on the Raychaudhuri equation, an approach described in detail in the next section (Sec. 0). This improves 
on previous explorations of the same issue which relied on a characterization of singularity avoidance based exclusively 
on the energy content of the universe through an extrapolation of the standard thermodynamics [20]. We thus provide 
an analysis of the implications of rainbow metrics on this issue which aims to account for all possible modifications. 

Then in Sec. |III| we examine the quantitative implications of the improved methodology of analysis here proposed for 
the cases of two particularly meaningful examples of rainbow metrics. This allows us to show that the improvements 
we advocate are not only valuable conceptually but also have important quantitative implications for the outcome of 
the analysis. 

The closing Sec. |IV| is devoted to a perspective on our proposals and results. There we also argue that our findings 
might suggest that the semi-classical approach which has been popular in the first explorations of rainbow cosmology, 
while being useful in getting hints about the early universe behavior, might not be appropriate all the way back to 
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the initial singularity, since in that regime quantum-gravity effects not captured by the semi-classical approach may 
be relevant. 

We use Planckian units such that H = Kb = c = 1. Notice that the deformations implied by MDRs and/or rainbow 
metrics do not weaken the role of speed-of-light scale c in physics, since that scale still is a relativistic invariant, the 
observer-independent value of the speed of massless particles in the infrared limit (the limit of small particle energy). 
For what concerns the speed of massless particles the only novelty one allows is that at high energies the speed might 
not be c. 


II. AVOIDING SINGULARITIES IN THE RAINBOW UNIVERSE 


The existence of singularities in GR is an intrinsic feature of the theory, as seen through the singularity theorems 
of Hawking and Penrose [21) , which are based on the notion of geodesic incompleteness. The analysis demands the 
use of the Raychaudhuri equation, i.e., the evolution equation for the expansion coefficient 9 = V of a congruence 
of curves defined by the velocity field V^. In the case of a congruence of geodesics it is given by 

0 = 2cu 2 - 2er 2 - l -9 2 - R llv V^V v , (1) 

where dot means derivative with respect to an affine parameter, 2w 2 = where lo^ v is the vorticity tensor and 

2er 2 = a^a^ 1 ' where is the shear tensor associated to the congruence. In the particular case of the Friedmann- 
Lemaitre-Robertson-Walker (FLRW) model with velocity field V^ = S(f, this equation reduces to the equation for the 
second derivative of the scale factor. 

The important point here is that if 9 does not change sign, then inevitably the congruence will shrink to a singularity 
if at some time to we have 9q < 0. This can be easily seen as follows: suppose that it is possible to construct 
as a gradient, that is, we define a congruence that is hypersurface orthogonal. In that case, the vorticity term w 2 
vanishes. The term related to the shear a 2 doesn’t change sign as it is a quantity orthogonal to V 11 . The last term 
is proportional to p + 3p if GR is assumed. Therefore, if the strong energy condition holds, that is p + 3p > 0, this 
term will never change sign as well. It follows that [2T 


dO 1 2 

— < — 9 2 => 9 
dr 3 


-l 


(t) >-t + 9 0 


( 2 ) 


Thus, if 0o < 0 for some value of r, then 9 1 —>• 0 and 9 —> —oo within the proper time interval r < 3/|0o|- 
In summary, upon assuming 


• the validity of GR and of the strong energy condition; 


• that gravity is strong enough to trap a region, so that one can have 9q < 0 at some 7o; 

• the existence of global hyperbolicity (which allows the definition of a global vector field with null vorticity); 

one finds that the space-time is geodesically incomplete, that is, the history of a particle may have an end or a 
beginning within a finite interval of proper time. These assumptions thus provide a set of conditions that imply the 
existence of a singularity. More generally, different combinations of the conditions above are possible. One could 
weaken one of them and assume stronger versions of the others to get different versions of the theorem [22] ■ If any of 
the conditions is not verified, then one can not conclude anything about the existence of the singularities and should 
look directly at the evolution equation, Eq. ([l]), and the invariants of the geometry constructed with it. Whereas 
previous studies of singularity avoidance in rainbow cosmology adopted simplified (and less conclusive) criteria, we 
shall here study geodesic (in)completeness through the Raychaudhuri equation 0 - 

In the context of a quantum-gravity theory, as mentioned, it is often useful to consider a semi-classical regime 
in which specific quantum-gravitational effects can be effectively described in terms of modified dispersion relations 
and the space-time can still be described as a smooth manifold, with non-trivial properties (low-energy relics of 
quantum-gravitational effects) encoded exclusively in an effective energy-dependent metric. 

A generic MDR can be cast in the form 

- f 2 (E)E 2 +g 2 (E)p 2 =m 2 , (3) 


1 We recall that a non-null shear can be a result of anisotropic expansion. Thus, as long as anisotropic spaces do not present vorticity, 
their singular behavior cannot be avoided based solely on the presence of shear. 
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where f(E) and g{E) are functions of the particle energy E and of the deformation scale Ep , which is expected to 
be of the order of the Planck energy. One must enforce of course that f(E) = 1 = g{E) in the limit E <C Ep, so that 
the standard special-relativistic dispersion relation is recovered in that limit 2 . 

From a MDR of the form (|3j) one can derive, in the case of a homogeneous and isotropic geometry described by a 
FLRW-like metric, the following line element [7] 

ds 2 = - ] dt 2 + a Sijdx'dx 3 . ( 4 ) 

f ( E ) 9 2 {E) 

Here a(t) is the standard expansion factor - the one describing the geometry probed by any particle with energy 
E <C Ep. When the particle energy is high enough that the modification of the dispersion relation, (J3|, cannot be 
neglected, the particle will ’feel’ the geometry described by Q, with an effective expansion factor a e // = a /g (see 
also discussion in j25j). 

The metric © is explicitly time-dependent because of a(t), but the time dependence comes also indirectly from 
the functions / and g. Indeed they depend on time through the energy, because of the redshift due to space-time 
curvature. Considering the normalized velocity field V M = f(E)5g we thus obtain the expansion coefficient 

where 6 e ff = 3a e f / /a e f / is defined in analogy with the FLRW case (where it reads 9 = 3 a/a) and whose evolution 
equation gives 3 


V/ 


0 2 eff 


= -y^(p + 3 p) - 7 


/ d eff 

f 3 ’ 


( 5 ) 


where dot from now on refers to derivative with respect to t. Before moving forward with the analysis it may be worth 
pausing here briefly to comment on our description of evolution in terms of the time variable t. At first sight one 
might be tempted to use a different time variable, say t ', obtained by absorbing in t the energy-dependent function /. 
However, evidently such a map from t to t! would not be a diffeomorphism, indeed because of the energy dependence 
involved 4 . Moreover, the new time variable t' would not correspond to the readout of any physical clock: the energy 
dependence postulates in rainbow-gravity proposals is meaningful exclusively when the relevant energy scale is the 
energy of a fundamental particle, but one cannot assign a similar notion of energy to a physical clock. It is also well 
understood IMETj that the relevant rainbow-gravity effects introduced for fundamental particles leave only negligible 
traces for macroscopic/composite bodies such as a physical clock. 

Note that if / is a well behaved function - it should be bounded for large values of energy in order not to give a 
degenerate metric and should smoothly tend to 1 in the low energy limit in order to recover the special relativistic 
dispersion relation - the behavior of 9 is reflected by the behavior of 9 e ff, and therefore we can explicitly account for 
the contribution of / to the evolution of the congruence by keeping it on the right-hand side of the above equation. 

The 0-0 component of the Einstein equations reads 


2 24tt 

e ff J2 P ’ 


and the continuity equation is given by 


p + 6 eff (p + p) = 0. 


( 6 ) 

( 7 ) 


2 Notice that this form of MDR is a deformation of the special relativistic one, not of the one one would have in a curved space-time. In 
this work we will assume that the generalization to curved space-time MDR’s can be made simply by taking into account redshift of 
energies, in the same spirit of [23]. A more rigorous treatment would require to consider deformations of a proper dispersion relation in 
curved space-time and the associated rainbow metric, but unfortunately this kind of studies is still at a very preliminary stage (see e.g. 

m\ 

3 This equation is not fully equivalent to the one presented by E3 and later by m because there is a missing term in these references 
proportional to f 

4 Note that even though we have left the standard setting of Riemannian geometry, an energy-dependent metric can be accommodated 
for example in a Finsler geometry context |181I19I where the invariance under usual diffeomorphisms is unaltered and energy-dependent 
coordinate transformations are not in general symmetries of the geometry. 
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Through Eq. ([5]) we see that there are basically two ways in which the initial singularity could be avoided: either 
the condition p + 3p > 0 is violated due to modifications of the equation of state p = ui(T)p - where the equation of 
state parameter becomes temperature-dependent - or the last term changes sign, which is possible depending on the 
form of the function /. 

The dependence on the temperature of the equation of state parameter arises from the modifications of the thermo¬ 
dynamics of the fluid being considered, as now the relation between energy and momentum is not trivial |28ti31j . The 
relevant thermodynamical quantities are computed using the density number of states with momenta in the interval 
between p and p + dp {521 which is defined by 

N(p)dp = ^ p 2 dp , (8) 

71 


and this relation gets modified when written in terms of energy due to the MDR, namely 


N(E)dE 


V 




E 2 dE, 


(9) 


where x' denotes derivative of x with respect to E. Therefore, one has a modification of the total average energy, 
defined for a gas of photons by 


U = 


E 


N(E) 


exp [E/T] - 1 


dE, 


( 10 ) 


and consequently of the energy density p = U/V. Other quantities like pressure also get modified, implying a deviation 
of the equation of state parameter u> = p/p in the UV regime. In the next section we study two interesting cases 
where these modifications are not negligible. 


III. QUANTITATIVE ANALYSIS OF SOME SINGULAR SOLUTIONS 

In the previous section we have discussed what are the relevant quantities to be studied in order to check for 
the presence/absence of singularities in a given rainbow universe. In particular it is possible to check directly the 
evolution equation governing the behavior of the effective scale factor a e ff which is given in ([5]) and the equation of 
state parameter should be computed explicitly once we have the rainbow functions. 

Before proceeding with our analysis, let us comment on some previous works dealing with the issue of singularities in 
the rainbow universe, since the way they treat thermodynamics is significantly different from what we do in this work, 
and leads to significantly different conclusions. The analysis presented in j20| has assumed the classical equipartition 
theorem to relate the average energy to the temperature without taking into account that a MDR should modify the 
thermodynamical properties of the fluid. That would only be true if / = g and for massless particles. Even if one 
wanted to claim that the corrections would be of second order in the Planck length ^T 2 , and therefore negligible 
[33] . the specific form of the MDR is needed to justify it. This means that the general criteria (i) and (ii) presented 
in [20} are not fully general, as the specific form of the MDR should be used from the beginning. This argument 
applies also to the fluid equation of state, which should also change due to the MDR. A proper analysis can be seen, 
for instance, in the context of black holes [Mj . 

Besides, in Refs. [20ll25ll33] the authors offer an average description of the metric Q, considered to be the effective 
metric probed by a radiation fluid with average energy E, but it should be appreciated that this leads to a picture in 
which the singularity is avoided on the average, while nothing can be said about probes with energies different from 
the average one. We show in what follows with explicit examples that this ambiguity is avoided when we consider a 
time dependence only through the actual cosmological expansion a(t), that is, the one defined for low values of energy, 
so that all the probes exhibit the same time behavior even though they perceive different geometries. 


A. First example 


In order to clarify the procedure we suggest to follow, and in order to emphasize the possibly different conclusions 
that we reach in comparison with previous analyses, let us now consider as an example the MDR that was studied in 
120) . where it was argued that it results in a non-singular solution. This particular form of MDR is also connected 
rather directly to a quantum-gravity picture, whose interest was first stressed by one of us in m ns]- It reads 


0 E/E P _ i 


E 2 = p 2 , 


E j Ep 


( 11 ) 
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so that the rainbow functions are given by 


e/e p _ 1 

f(E) = E/Ep , 9(E) = 1. (12) 

Firstly we study the thermodynamic properties of a gas of photons characterized by such a dispersion relation, then 
we show how the modifications in the statistical quantities should be considered in order to obtain explicitly the 
evolution of the scale factor, which turns out to be singular contrary to what was claimed in [201 . 


1. Modified thermodynamics 


For an ensemble of particles, the density number of states with momenta in the interval between p 
given by Eq. ([8]). Therefore, it can be written in terms of energy by using the MDR governed by (12) 


and p + dp is 
as: 


N(E)dE = ^-E 2 P (e E/Ep - 1 ) 2 e E/Ep dE. 

7 T Z 


(13) 


The equation of state parameter is defined in terms of the average energy density and the pressure according to 


p _ / In [1 - e~ E / T }N(E)dE 

P / exp [E/T]- 1 E(E)dE 


(14) 


The above expression can be numerically integrated to give us the equation of state parameter in terms of temper¬ 
ature, as shown in Fig jlj We see that u> decreases from its usual value 1/3 as the temperature increases and it quickly 



FIG. 1: Plot of the equation of state parameter as a function of temperature (Tpi is the Planck temperature). 

approaches a constant asymptotic value, namely, u> —> 0 in the UV regime. It is interesting to note that in this regime 
the radiation fluid is effectively behaving as dust. 


2. Cosmological evolution and existence of singularity 

In this example, g(E) = 1 and thus a e ff coincides with the actual scale factor. Besides, according to what has 
just been shown, we can assume that a; ~ 0 in the UV regime. Then, the continuity equation 0 gives us the energy 
density as a function of the scale factor as p oc a~ 3 . Besides, assuming that the time dependence of physical energy 
and momentum scales with the scale factor as p = k/a, k being the comoving wavenumber, we can obtain from © 
the relation 


k 

Epa 



E = Ep In 


(15) 
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so that the rainbow function f(E) can be written as a function of the wavenumber k (associated to energy E) and of 
the scale factor: 


f(a) = 


In 


k/ ( aEp ) 

(e^ + 1 ) 


Plugging p(a) and /(a) into the first-order Friedmann equation, Eq. § , we obtain 

8-7T p p Ep / - ] ( k 


a = 


\fa\x 


\Epa 


+ 1 


(16) 


(17) 


where p p is a constant of integration coming from the continuity equation. Note that this equation is valid only in 
the UV, for which the equation of state parameter takes the value u> = 0. 

This equation leads to a singular solution as we can explicitly see in Fig. [2] where we plot the solution of this equation 
for different values of the wavenumber in comparison with the standard radiation-dominated FLRW solution. 



t 

p 


FIG. 2: Comparison between different solutions for the scale factor where subscript p refers to Planckian 
is plotted the standard radiation-dominated FLRW solution; in green, dark blue and red are, respectively, 
modified Eq. (171 for k = l/107? p , k = 3/10 Ep and k = Ep. Note how the modified evolution slows down 
singularity when compared to the standard FLRW case. Here t p = l/^/87rp p /3. 


quantities. In black 
the solutions to the 
as it approaches the 


We see that the modified evolution slows down near the singularity such that the higher the energy the slower 
the approach. Indeed, the concavity of the curve changes for large values of energy indicating that quantum gravity 
effects are somehow smoothing the singularity. However, this is not enough to avoid the singularity, although it 
might be suggestive of the fact that a full quantum-gravity picture, admitting our effective-theory description as an 
approximate description, could actually avoid the singularity. 

Another way of seeing this phenomenon of slowing-down is through the comparison between the Hubble parameter 
coming from Eq. (17]) (defined as Hk = a/a) and the corresponding quantity in the FLRW radiation-dominated 
universe as the scale factor approaches zero. We remark that the Hubble parameter is invariant under spatial re¬ 
scaling transformations and therefore the slowing-down of the modified solutions is independent of the specific values 
we set for the comoving wavenumber. Thus, we have 

—r - oc In ( —p —b 1^ => lim —p = 0. (18) 

H \E p a J a —>o H 

That is, even though both quantities diverge as a — > 0, their ratio goes to zero for fixed values of fc, meaning that the 
expansion rate for the FLRW case diverges faster than in the modified picture. 

Note that in principle the singularity could be avoided from the arguments of section [h] The term /// in Eq. ([ 5 ]) 
can be explicitly written as 

/ = _ (_?L _ I\ ( 1 \ Oefl 

f \e x — 1 x) \l + Epa/k) 3 ’ 


( 19 ) 
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where x = E/Ep. The term in the first brackets is always bigger than zero for all values of x. Plugging this 
expression back into Eq. ([5| it is possible to see that it contributes with an opposite sign and thus, in principle, it 
could prevent the divergence of d e ff- The fact that our semi-classical description, based on an “effective” rainbow 
metric, is unreliable as the scale factor goes to zero is seen through the fact that arbitrarily high values of energy can 
be achieved and therefore the rainbow function and consequently 0 e // diverge. 5 


B. Second example 

Let us now consider another much studied example of the MDR for massless particles, given by the relation 

E 2 =p 2 [l + (Xp)% (20) 

so that the rainbow functions are / = 1 and g(p) = [1 + (Ap) 4 ] 1 / 2 . A is a deformation parameter with dimensions of 
inverse of energy. This MDR is found to provide an invariant power spectrum for cosmological perturbations and also 
gives the expected spectral dimension reduction found in different approaches to quantum gravity [36] : the UV value 
of the spectral dimension is two. Again, we first analyze its thermodynamic properties and then we study this fluid 
in the modified FLRW geometry showing that the cosmological evolution also contains a singularity. 


1. Modified thermodynamics 


The density of states is given as before by Eq. (Jsj) where we now keep the dependence on p. As in the first example, 
the equation of state parameter varies, departing from the usual value 1/3 and approaching to = 1 in the UV limit, 
as seen from the numerical integration of the expression 


p f In [1 — e E{ ' P ^ T ]N(p)dp 

9 f exp [E(r)/T]-l N (P) dp 

in Fig. [| 


( 21 ) 





FIG. 3: Logarithmic plot of the equation of state parameter as a function of T/Tp with A = 10 5 Lp, with Lp the Planck length. 

It is intriguing to observe that the UV value of the equation of state parameter matches the one of a radiation 
fluid in two spacetime dimensions, which is also the value of the UV spectral dimension that can be deduced from the 


5 It could be the case that a MDR is capable of providing a bouncing solution such that the bounce would occur at an energy scale 
several orders of magnitude smaller than the Planck energy. Then an analysis based on the full quantum gravitational theory wouldn’t 
be strictly necessary. 







dispersion relation (201 


Moreover, also the energy density behaves in the UV as if the radiation fluid was two- 


dimensional. To show this we perform a numerical integration and find the total energy as a function of temperature, 
through the expression 


U = 



N(p) 

exp [E(p)/T\ - 1 


dp, 


which yields the curve shown in Fig. [4j 


( 22 ) 



FIG. 4: Average energy as a function of T with A = 10 5 Lp. 

It is possible to see that the system seems to behave as a conventional one (usual dispersion relation), but in two 
space-time dimensions. By fitting the curve one can easily check that in the UV it goes like T 2 . If we assume that 
the total energy of the system goes with U oc T d+1 , where d is the number of spatial dimensions [2H[ [22], we see that 
our system is effectively behaving like one with the usual dispersion relation in one spatial dimension. 


2. Cosmological evolution and existence of singularity 

Again, the time dependence of the function g comes only from the time dependence of the actual scale factor 
a(t), thus it is the same for all particles with different momenta. This means that even if the probes are capable 
of perceiving only the effective scale factor a e //, and the latter is dependent on the wavenumber of the probe, the 
presence of a singularity will be a general feature of the space-time. 

In the case we are considering / = 1, the only possibility to avoid a singularity according to Eq. © would be a 
non-conventional behavior of the fluid such that it could violate the energy condition p + 3p > 0. As we saw in the 
last subsection, this is not the case for the MDR considered, as the equation of state parameter asymptotically goes to 
the constant value u> ~ 1, resulting in a time dependence of the form a e ff(t ) oc t 1 / 3 . Therefore, this rainbow universe 
is singular, as can be explicitly seen in Fig. [5] 

We see that the modified evolution seems to be slower when compared to the radiation dominated FLRW case for 
most of the evolution of the universe, which is an indication that it could bypass the singular behavior. However, as 
we follow it backwards in time it actually reaches the singularity even more abruptly than in the standard case. Again, 
the analysis indicates that the full quantum gravity should be applied in the far UV regime, where the semi-classical 
picture is no longer an accurate description of the quantum universe. 


IV. CONCLUSIONS 

As some previous studies we here worked within the legitimate assumption (legitimate but only one of many possible 
hypotheses) that a rainbow metric Q could be used for a reliable approximate effective-theory description of some 
quantum-gravity effects. Our strategy of analysis provides an improvement on previous analogous studies in several 
ways. We have shown (and used the fact) that a proper analysis of the possible avoidance of the initial singularity 
should take into account the modified thermodynamics implied by the MDR, as it can drastically change the time 
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t 
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FIG. 5: Comparison between the standard evolution of the scale factor in a radiation dominated FLRW 
which g(E) = 1 and the evolution of the effective scale factor in the case of the MDR given in Eq. (201 
value of the wavenumber. 


universe (in red) for 
(in blue) for a fixed 


behavior of the scale factor. Our approach is also not based on an average description as the one introduced in [551 
and later used in 150], so it is more conclusive in that our findings apply also to probes with energy different from 
the average one. We have argued that it is possible to provide a consistent analysis of the issue considering a time 
dependence that comes only from the cosmological expansion. We have provided explicit computations showing, in two 
different pictures, that our concerns are not merely conceptual but rather have significant quantitative implications 
(to the point that pictures found to avoid the singularity in other more simplified approaches are not found to avoid 
the singularity within our approach). Still, from the physics perspective the most exciting feature is the one we found 
in the first of the two specific models here analyzed, where the approach to the singularity is slowed down by the 
MDR. As also supported by some of the points made here, such cases where singularity approach is slowed down 
within the effective-theory description based on a rainbow metric may well ultimately be part of a quantum-gravity 
picture in which the singularity is avoided, even though extrapolating the rainbow-metric approach to very small scale 
factors we do still find the singularity. This is plausible because the onset of features slowing down the singularity 
approach occurs in a regime which may well be within the reach of the effective-theory description based on rainbow 
metrics, whereas at even earlier times (smaller values of the scale factor) the reliability of this effective description 
most evidently breaks down. The additional quantum-gravity effects that one should take into account at that stage, 
intervening in a scenario where singularity approach has already been slowed down, might well provide the needed 
additional structures for avoiding the singularity. 
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